Spin exchange interaction with tunable range between graphene quantum dots 
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We study the spin exchange between two electrons localized in separate quantum dots in graphene. 
The electronic states in the conduction band are coupled indirectly by tunneling to a common 
continuum of delocalized states in the valence band. As a model, we use a two-impurity Anderson 
Hamiltonian which we subsequently transform into an effective spin Hamiltonian by way of a two- 
stage Schrieffer- Wolff transformation. We then compare our result to that from a Coqblin-Schrieffer 
approach as well as to fourth order perturbation theory. 
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I. INTRODUCTION 

Spins in quantum dots (QDs) are under intense in- 
vestigation as a possible realization of a quantum bits 
(qubits) ji Among the currently most advanced solid- 
state structures are top-gate patterned two-dimensional 
electron gases in GaAs heterostructures.— However, in 
this host material hyperfine interaction between the spin 
of the electron and that of the atomic nuclei of the host 
material leads to relatively short coherence times. A 
promising way to circumvent this problem is the use of 
carbon as a host material for spin qubits. Natural car- 
bon comprises 99% of the carbon isotope 12 C which has 
no nuclear spin. This gives carbon based devices the 
advantage that decoherence due to hyperfine interaction 
is suppressed by the small abundance of nuclear spins. 
Carbon is also a relatively light element, therefore spin- 
orbit coupling is expected to be weaker than in GaAs. 
One can expect a significant improvement of spin co- 
herence times in carbon based structures. Graphene, a 
truly two-dimensional carbon-based crystal 3 - - — seems to 
be an ideal host material for spin qubits^. It naturally 
created a perfect confinement of electrons in one dimen- 
sion. Moreover, in contrast to carbon nanotubes ; 10 ' 11 the 
ability to lithographically pattern graphene allows for a 
deterministic device preparation, which is necessary for 
scalability^ Graphene has a very interesting electronic 
structure with a gapless and linear dispersion around the 
Fermi energy. Furthermore, the electronic eigenstates 
carry an additional internal degree of freedom, dubbed 
pscudospin, which is always aligned with the direction 
of the momentum. These properties imitate the be- 
havior of rclativistic chiral massless Dirac particles^^ 
These rclativistic-likc properties lead to the phenomenon 
of Klein tunneling ) 12 : 13 which actually prohibits any elec- 
trostatical confinement of electrons, i.e. prohibits the 
formation of quantum dots. 

Among the most promising ideas to overcome Klein 
tunneling is to use graphene nanoribbons or constrictions 
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FIG. 1: One-dimensional quantum dot array on an armchair 
graphene nanoribbon (drawing not to scale). Due to the rib- 
bon structure, the dispersion relation of graphene can exhibit 
a gap, which scales inversely with the ribbon width. This gap 
allows for electrostatic confinement of electrons in quantum 
dots. As the bandgap is small compared to regular semicon- 
ductors, the spin exchange mechanism between the quantum 
dots is not dominated by RKKY-type processes alone, and 
superexchange processes contribute significantly. 



instead of extended graphene as host material ^ 14 : 15 see 
Fig. [TJ In clean graphene nanoribbons with armchair 
boundaries, the additional confinement can lead to the 
opening of a small energy gap at the Fermi energy*^ The 
size of the gap is indirectly proportional to the width of 
the ribbon. In the presence of such a gap, the pseudo- 
rclativistic behavior of the charge carriers is lost, and 
the material resembles a regular gaped semiconductor, 
enabling electrostatic confinement £ Interestingly, exper- 
iments observe the formation of a gap irrespective of the 
boundary condition ^2. As a side effect, the sharp edges 
also lift the valley degeneracy in bulk graphene, which 
could suppress the Hcisenberg spin interaction between 
the quantum dots4 

Following Ref. 0, we consider a system as shown in 
the upper part of Fig. [T] where several electric gates are 
placed on top of an armchair nanoribbon. By applying a 
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gate voltage, the dispersion relation of the material below 
the gates can be shifted in energy, see the lower part 
of Fig. [TJ If at a certain energy extended states could 
exist in one section, but not in the neighboring ribbon 
sections, additional size quantization along the ribbon 
leads to single localized states. In the energy interval 
above and below the gap, the states are extended and 
form a continuum. With a suitable adjustment of the 
chemical potential, the localized states can be filled with 
one electron each, forming a one-dimensional array of 
qubits. In the following, we calculate the Heisenberg spin 
interaction between two such spin qubits. 

The RKKY-interactior>i£ between localized magnetic 
moments is well discussed for extended graphene, see 
Refs. GUH and after the ex peri mental discovery of 
graphene revisited in Refs. I21K251 Here, we study a 
graphene nanoribbon where a gap opens at the Fermi en- 
ergy. Due to this gap, the spin exchange problem we are 
interested in resembles less the one in extended graphene 
but more the case in ordinary semiconductors in reduced 
dimensionalities,—"— like quantum wells^ 9 - or quantum 
wires j 3 - - However, compared to conventional semiconduc- 
tors, the band gap can be unusually small, on the order 
of 1 — 100 meV. Therefore, by applying gate voltages, 
one can realize an arbitrary alignment of the quantum 
dot energy level relative to the valence and conduction 
band. Due to the proximity of the band edges of valence 
and/or conduction band, the band edges must be taken 
into account^ 



turn dots are occupied with one electron each to form a 
spin 1/2 qubit. 

The second part of Ho models the contacting contin- 
uum of states as a large reservoir of noninteracting elec- 
trons. Here Cka, c\. a denote the annihilation and creation 
operators for electrons in the continuum with (orbital) 
quantum numbers k and spin a. The continuum is as- 
sumed to be unpolarized and at zero temperature (filled 
valence band). The tunneling Hamiltonian Ht describes 
spin-conserving tunneling between the two dots and the 
continuum. The tunnel amplitudes Uk(r.i) depend on the 
dot and the continuum quantum numbers i nd k, as well 
as on the position r.; of the quantum dots. However, 
the exact form of tik{ri) depends on the system under 
consideration. 

The two-impurity Anderson model and the strongly 
related Anderson lattice model is extensively dis- 
cussed in the literature in the context of rare earth 
compounds^"— dilute magnetic semiconductors^ high 
temperature superconductors.—"— Numerous theoreti- 
cal techniques were used to study these models, in- 
cluding variational wave functions ; 45 ' 46 equations of 
motion ) 36 ! 37 finit o 34 ' 35 i 42 i 47 ' 48 and infinite perturba- 
tion expansions^ 9 - Bethe-ansatz studies^ higher-order 
Schrieffer- Wolff transformations ; 43 ' 44 ' 51 " — or, most com- 
mon, a mixture of a first-order Schrieffer- Wolff trans- 
formation and perturbation theory ; 33 ' 40 ' 41 ' 56 " — originally 
proposed by Coqblin and Schrieffer in Ref. [H. 



II. TRANSFORMATION OF THE 
HAMILTONIAN 

A. Model 

We model the system as two Anderson impurities, 
which arc both in contact with a common energy band. 
The Hamiltonian for our model is 

H = Ho+Ht, (1) 

i=l,2 \<r=f4 / 

+ Y. £kC l^ka, (2) 

ka 

H T = J2 (^(^) c L a i<r + 4(^)«t c ^) • ( 3 ) 

ika 

The first part of Hq describes the two independent quan- 
tum dots i = 1,2, containing one electronic level each. 
The fermionic operators a ia and a\ a create and annihilate 
electrons on dot i with spin a. Due to the low electro- 
static capacity of a quantum dot, double occupation of 
one individual dot is associated with a substantial charg- 
ing energy EZj. For quantum computation applications 
one is interested in a parameter regime where the quan- 



B. Schrieffer-Wolff Transformation 

Following Refs. l5lM54l we use a two-stage or nested 
Schrieffer-Wolff transformation to derive an effective spin 
Hamiltonian. In contrast to previous works we do not as- 
sume equal energy levels in the two Anderson impurities, 
as the confinement of the quantum dots can be modi- 
fied individually. By keeping track of the dot indices it 
is also possible to identify different physical processes in 
the final result, and enables us to compare our result to 
higher-order perturbation theory. 

The Schrieffer-Wolff transformatio n 38 ' 59 ' 62 " — is based 
on a canonical transformation of the Hamiltonian, 
H (i) = e iS H e -iS = H+[iS,H] + ±[iS, [iS,H]}+.... The 
division of the Hamiltonian H into a free Hamiltonian Hq 
and a small perturbation Ht, allows us to choose a trans- 
formation Si, fullfilling the relation [iSi,Ho] = —Ht 
and leading us to the effective Hamiltonian H^ = Ho + 
\[iSi, H T ] + \[iSi, [iSi ,H T }] + ± [iSi , [iSi , [iSi , H T ]]] + ... 
where the lowest-order tunneling term is canceled (note 
that the coefficient has the general form 1/nl — l/(n + 
1)! = l/[(n - l)!(n + 1)!]). Since S also has to be of 
first order in the tunneling amplitudes, Si on Ht, the 
interaction now appears (at least) in second order. 

By a subsequent Schrieffer-Wolff transformation with 
the generator S2 fulfilling [iS2,Ho] = — h [iSi , Ht] , also 
the second order interaction term can be removed. Note 
that now, S2 oc H T . At the end we project the result- 
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ing Hamiltonian on the subspace where both quantum 
dots are occupied by one electron. As all odd-order in- 
teractions do not conserve the occupation numbers of the 
quantum dot, they can be neglected, as they will be pro- 
jected out at the end of the calculation. Combining both 
steps, we arrive at the effective Hamiltonian 

H {2) = H + hiS2, [iSx,H T ]] + l[iS 1: [iS u [iSi,H T ]]], (4) 

4 o 

where corrections in sixth and higher orders in the tun- 
neling amplitudes have been neglected. After project- 



ing out the continuum degrees of freedom, and in ad- 
dition to unimportant level renormalizations, which we 
do not discuss, we find a Heisenbcrg-like interaction 
J S± ■ S2, which couples the two quantum dot spins 
Si = J2a8 a ta (T °>P a i0i consistent with Refs. l5ll|52l A 
detailed calculation is presented in the Appendix. Af- 
ter non-trivial regrouping of terms^, one can separate 
the spin-interaction into parts originating from different 
virtual tunneling processes defined by their intermediate 
virtual quantum state with the explicite shape: 



J 

Ji 
■h 
■h 

Ji 



k.q 

\ rift — n q ( 1 



Sk — £l e q — £i—UiJ Eft — Eg \Sq — £2 £k — £2 — U; 



1 



1 



1 — n„ 



1 



1 



£k — £l £q — £2 J £k + £q — £l — £2 \£q — £\ £k ~ £2 
1 1 \ +71, 



1 



1 



£l + Ui-£ k £2 + U 2 - £ q J £l + Ui + £ 2 + U 2 - £k - £q \£l + Ui - £ q £2 + U 2 -£kJ 7 

1 1 \ -n n ( 1 1 \ 1 +1 1 



£k — £\ £k — £2 — U2 J £2 + U2 — £1 \£ q — £1 £ q — £2 — U2 J £k — £\ £2 + U2 — £1 £q — £1 



(5) 

(6) 
(7) 
(8) 

+ (lo2). (9) 



The first term Ji resembles an RKKY-interactionJ^ 
The interaction is mediated by a virtual particle-hole ex- 
citation in the electron gas, see Fig. EJa). Therefore the 
energy of the intermediate excitation is given by £ q — 
The second and third contributiosn to the spin-spin in- 
teraction originate from virtual two-particlc(hole) excita- 
tions in the continuum Fermi sea, see Fig. Et^b) and (c). 
Two electrons tunnel coherently from or to the quantum 
dots. Thus, as intermediate virtual states, the two quan- 
tum dots are both doubly occupied (empty). Afterwards, 
the electrons tunneling crosswise back, interchanging the 
spins of the quantum dots. This process leads to the 
interactions J 2 and J3. 

Finally, the last contribution J4 is caused by the pos- 
sibility of a direct tunneling of one dot electron to the 
other dot. The virtual intermediate state is therefore 
one double occupied dot and one empty dot. The tun- 
neling can happen through filled as well as empty states 
in the electron gas, see Fig. ^d) 



C. Relation to Coqblin-Schrieffer Model 

In Ref.dHCoqblin and Schrieffer presented their widely 
used approac h 33 ' 50 i ' 66 to the two-impurity Anderson 
model. They performed a single Schrieffer- Wolff trans- 
formation, and project the resulting Hamiltonian of the 
single-occupied impurities. With this, they transform the 



two individual Anderson models into two individual s-d 
models or Kondo impurities. The spin of one quantum 
dot Si is coupled to the bath electron spins by 

-ffKondo = ^2 JkqSi • C^O •„„, C q(T , , (10) 
kq,aa' 

see Appendix A for details. By treating that Hamilto- 
nian in second order perturbation theory, they compute 
a RKKY-like spin-spin interaction^ of the form 

£«ft^^-s 2 . (ii) 

kq £k "1 

Even though Eq. (TTTj) captures the basic features of 
Eq. ([6]), it is an expansion inconsistent in the order of 
the tunnel amplitudes. First, the initial Schrieffer- Wolff 
transformation generates not only terms of second order, 
but also the term [iSi, [iSi, [iSi, Ht] which contributes in 
fourth order,— see Eq. ([4|. Actually contributions from 
this higher-order term cancel several parts in Eq. (jTTJ) , 
leading to Eq. ©. By truncating the transformation at 
second order, these contributions to the spin interaction 
are lost. 

Second, the initial Schrieffer- Wolff transformation does 
not only generate the Kondo-Hamiltonian, but also terms 
in second order of the form a\ 02, which describe direct 
tunneling between the two quantum dots^ In a subse- 
quent second order perturbation theory, these terms also 
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(c) Two-continuum-particle excitation 
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(d) Direct tunneling 




FIG. 2: Several virtual tunnel processes contribute to the 
spin-spin interaction between the dots. These processes can 
be classified by the intermediate state of the system. While 
particle-hole excitation (a) leads to an RKKY-like interac- 
tion, the processes (b)-(d) are usually summarized as superex- 
change. 



lead to a inter-dot spin interaction. In Ref. [H these in- 
teractions are neglected due to the premature projection 
of the result of the Schricffcr- Wolff transformation on the 
single occupied dot subspacc. 

Interestingly, in the limit of energy levels far away 
from the Fermi energy, i.e., when one assumes that the 
spin coupling J % kq approaches a constant, the Coqblin- 
Schricffcr approach generates the correct result. How- 
ever, nowadays the Anderson model is extensively used 
to describe quantum dots . 30 i 57 i 61 In contrast to rare earth 
compounds or true atomic systems, the typical energy 
scale of the quantum dot level spectrum is several or- 
ders of magnitude smaller. Therefore in these artificial 
systems the application of the Coqblin-Schricffcr model 
needs to be handled with care. 



D. Relation to fourth order perturbation theory 

Starting from the two-impurity Anderson model, one 
can also derive a fourth order dot spin-spin interaction 
by perturbation theory , 34 ' 35 ! 40 ' 47 ' 48 ' 67 with or without di- 
agrammatical help. The perturbation approach nearly 
reproduces our results Eq. (jB])-© with one exception: 
the structure of the Fermi functions. Via perturbation 
theory one would expect for example that the two-lead- 
particle excitation, see Fig.f2Jc), only happens if the two 
electron gas states k and q are empty, therefore the spin 
coupling J2 should be proportional to (1 — nfc)(l — n q ). In 
contrast, the contribution from a Schrieffer- Wolff trans- 
formation is proportional to (1 — n q ), and independent 
of rife. By counting the operator commutators, one can 
directly determine that the spin-spin coupling derived by 
a two-stage Schrieffer- Wolff transformation can not gen- 
erate terms, which contain four lead operators, which 
would be necessary for a product term like n^riq. The 
reason for this discrepancy between fourth order pertur- 
bation theory and Schrieffer- Wolff transformation lies in 
the procedure of integrating out the lead degrees of free- 
dom, i.e., by the replacement of the thermal average of 
lead operators (c^c/Jth by the Fermi function nu- In the 
case of the perturbation theory, the operators c\ , cj, re- 
fer to the bare unperturbed electronic states of the lead, 
i.e., one assumes, that the lead is actually not perturbed 
by tunneling. After the Schrieffer- Wolff transformation, 
the lead operators refer to new lead states, which are 
hybridized with the localized dot states. By performing 
the thermal average, one therefore assumes that these 
new hybridized lead states are in thermal equilibrium, 
not the bare lead states. Therefore, it is not surprising 
that the results of a Schrieffer- Wolff transformation and 
perturbation theory differ. However, it is surprising that 
one can express the result of the Schrieffer- Wolff trans- 
formation in the same functional form one would expect 
from perturbation theory, except for the Fermi functions. 
Only due to this structure of terms one actually can, in 
the spirit of Feynman diagrams, assign virtual processes 
as shown in Fig. [2J For this reason, the grouping of terms 
in Eq. ([6])- ([9]) is physically plausible but to some extent 
arbitrary. It has already been observed by Rudcrman and 
Kittel^ that by assuming certain symmetries of the tun- 
nel amplitudes, the Pauli exclusion principle can actually 
become in part unimportant. 



III. APPLICATION TO GRAPHENE 
NANORIBBON QUANTUM DOTS 

Up to now, the computed result in Eq. is gen- 

eral for the spin coupling of two qubits by a common 
continuum of states labeled by the indeces k and q. In 
the following, we specify this continuum to the electronic 
structure of a graphene nanoribbon aligned along the y- 
dircction. 
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A. Band structure 

Bulk grapheme has two independent Fermi points at 
the momenta K and K' in reciprocal space, generat- 
ing the valley degeneracy. Due to the armchair bound- 
ary conditions of the ribbon, the propagating wave states 
with momentum K + k and K' + k are coupledi 16 i 68 The 
confinement in x-direction leads to a further quantization 
of the transverse wave vector k x = k n = (n ± l/3)n/W 
with W denoting the ribbon's width and neff. There- 
fore the continuum states can be characterized by the 
subband index n and the momentum component k y = k 
along the ribbon. Close to the Fermi energy, the disper- 
sion relations becomes 

£k,n = ftvp^k 2 + fcf, (12) 

with the Fermi velocity vf of graphene. This dispersion 
resembles the dispersion of a massive rclativistic particle. 

The transverse confinement determines the energy gap 
2s g = 2hvpko, which scales inversely with the rib- 
bon width. 17 Due to this gap, electrons can be con- 
fined by electrostatic gates, in analogy to conventional 
semiconductors 4 We assume the applied electric poten- 
tial to be independent of the x-coordinate (see Fig. [!} , 
therefore the band index n is conserved. Therfore we only 
need to consider the continuum subband with the same 
band index as the bound state(s). Even if this symmetry 
is broken, the generalization to multi-subbands is straight 
forward*^ As a further simplification, we assume that by 
applying electrostatic gates, the dispersion relation of the 
extended states is still described by Eq. (fl"2j). 
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FIG. 3: If the quantum dot level lies close to the valence band, 
and the charging energy is smaller than the band gap, direct 
tunneling processes will dominate the spin exchange between 
the quantum dots. 

the overlap of the wave functions^^ leads to a tunnel 
amplitude of the form 

p —iK-ri _i_ p—iK'-n 
Kkirx) = he- lk ^ ±= . (13) 

The spin coupling therefore will always contain a con- 
tribution which oscillates on inter-atomic distances, and 
one contribution, which varies on the length scale of the 
envelope wave function. As the quantum dots under con- 
sideration are not spatially defined with lattice site pre- 
cision, we expect that the oscillating contribution to the 
spin exchange will average out. 

C. Spin-exchange range 



B. Tunnel amplitudes 

The spin exchange is proportional to the product of the 
four tunnel amplitudes t\ k {r\)t 2 k( r 2)ti,q( r i)t2, q ( r 2) ■ in 
analogy to most cases studied in literature, we assume 
that the amplitude of the overlap of the bound states 
and the extended states does not explicitely depend on 
the momentum k of the extended states. This assump- 
tion is valid if the wave function of the bound state is 
localized on a length scale smaller than the wave length 
of the extended state. In this case, one can approximate 
the localized wave function as a delta-function. How- 
ever, in quantum dots in semiconductors in general, and 
in particular in the vicinity of a band edge, this assump- 
tion may not be valid. As fc-independent tunnel ampli- 
tudes lead to a shorter spin exchange range, the spin 
exchange range derived in the following can be seen only 
as lower bound. Although the magnitude of the tunnel 
amplitude does not depend on k within this approxima- 
tion, the fact that the two quantum dots are separated 
in space gives rise to a relative phase. While in ordinary 
isotropic metals this phase is simply given by e Vi , in 
graphene the valley degeneracy has to be taken into ac- 
count. In nanoribbons, the energy eigenstates are phase- 
locked superposition of states of both valleys. Therefore, 



Which of the virtual tunnel process (see Fig. [2]) dom- 
inates the spin exchange between two quantum dots de- 
pends on the alignment of the dot energy levels, band 
gap, and edges. Roughly speaking, the virtual process 
requiring the lowest excitation energy will dominate. As 
we assume the graphene nanoribbon to be nearly un- 
doped, the Fermi energy of the system lies within the 
band gap. Therefore the valence band is entirely filled, 
and the conduction band is empty. The RKKY-likc ex- 
change interaction via a particle-hole excitation in the 
continuum will be suppressed by the band gap 2e g , and 
superexchange processes will dominate. In the following, 
we will discuss the scenario shown in Fig. [31 where the 
quantum dot level lies close to the valence band, and the 
charging energy is smaller than the band gap. 

In this level alignment, direct tunneling processes via 
the valence band dominate the spin exchange. The re- 
sulting integrals can be computed by Cauchy's integral 
formula. The exchange due to direct tunneling turns out 
to be 

M 2 M 2 1 (S2 + U2) 2 

AAE 2 e 2 + U 2 - si e 2 g - (e 2 - U 2 ) 2 
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with the quantum dot distance Ar = \n — ra|. The 
energy AE = hvp/L is the energy splitting of the con- 
tinuum states, with L being the length of the ribbon. As 
the tunnel amplitudes U decrease with the real space par- 
ticle density of the band states with 1 / \^X, the strength 
of the spin exchange is independent of the overall length 
of the ribbon. 

The strength of the spin coupling driven by direct 
tunneling diverges within 4 th order, if the single occu- 
pied state of one dot becomes resonant with the dou- 
ble occupied state of the other dot. The range A = 

hv F /2^Je 2 g - (e~ 2 + U 2 ) 2 of the coupling Eq. flU]) on the 

other side is controlled by the energy separation between 
the double occupied state of the one dot, and the valence 
band edge at energy — s g . If one assumes, that the single 
occupied states of the two dots are close to the band edge 
(s g w £;), the exchange range scales as Hvf/ yJ$s g Ui. For 
a graphene nanoribbon with width of 50 nm and a quan- 
tum dot with a charging energy of 4 meV^ this length 
is of order 50 nm, i.e. comparable to the quantum dot 
length. 

If one assumes that, in analogy to a capacitor, the 
charging energy of a quantum dot scales inversely with 
its area^ then the range of the spin exchange interaction 
scales linearly with the width of the nanoribbon. 



D. Further considerations 



IV. CONCLUSIONS 

In this paper, we discussed the spin exchange be- 
tween localized states which are only indirectly coupled 
via a continuum of states. Using a 2-stage Schricffcr- 
Wolff transformation, we transformed a 2-impurity An- 
derson model into an effective spin Hamiltonian. Based 
on our result, we discussed the validity of the Coqblin- 
Schrieffer approach to this problem. Furthermore, by a 
re-ordering the terms, we were able to directly compare 
the Schrieffer- Wolff result to perturbation theory, and 
observe distinct differences that originate from different 
assumptions on the continuum of states. 

As an application of the formalism developed here, we 
discussed the spin exchange interaction between electro- 
statically confined quantum dots in a graphene nanorib- 
bon, as shown in Eq. (|14p. As a lower bound, we derive 
a range of this spin exchange of the order of the nanorib- 
bon width. However, the dot energies can be adjusted in 
such a way as to extend the exchange coupling to longer 
distances. 
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The lower bound for the spin exchange range between 
quantum dots in a graphene nanoribbon is given by the 
ribbon width. For this result we considered only virtual 
tunnel processes via free continuum states. In addition, 
also direct tunneling between the quantum dot states can 
occur. If the quantum dot level approaches the band gap, 
the bound electron leaks further and further into the bar- 
rier due to the weakening of the confinement. Therefore 
it can happen that two neighboring quantum dot wave 
functions can acquire a non- vanishing overlap, and direct 
tunneling becomes possible. Direct tunneling is accom- 
panied also by a spin-exchange interaction^ 

Furthermore, we assumed so far that the graphene 
nanoribbon is infinitely long. This assumption is hid- 
den in the approach to treat the continuum states with 
momentum k independent on the state —k. However, if 
a finite ribbon length leads to a defined phase relation 
of the forward and backward propagating states, then 
one part of the tunnel amplitude will become entirely 
independent on the momentum k, and therefore on the 
distance Ar& (cf. also the discussion of section UlIBI ) 
If this case, the range A of the spin exchange is not de- 
termined by the dephasing of the exchange contributions 
of different states within the Fermi sea, but only on the 
phase coherence length of the extended states. 



Appendix A: First Schrieffer- Wolff Transformation 

The generator of the first Schrieffer- Wolff transforma- 
tion must fullfill the relation [iS\,H^\ = —Ht- As the 
part of the Hamiltonian Hq only contain number opera- 
tors, i.e., is quadratic or quartic in the fcrmion operators, 
one can deduce, that iS must have the same structure as 
the tunnel Hamiltonian Ht, up to a prcfactor contain- 
ing either constants or number operators. With such an 
ansatz, the generator of the first Schrieffer- Wolff trans- 
formation is found to be 



iSi = 

E 



(Al) 



ti,k 



£k - £i - Uillig 
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1 - Uja 



s k - Si - Uin lS 



h.C. 



This transformation removes the interaction term of first 
order in the tunneling amplitude, but instead generates 
higher order interactions starting in second order of t. 
The new interaction Hamiltonian H T ^ = j[iS\, Ht] + 
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0(t 3 ) becomes 



Appendix B: Second Schrieffer- Wolff Transformation 



Bt = E J ^ ■ s ^ (A2) 

i.kq.a 

+ ~ E ^i3 a lcr a jc 

+ dot double empty /filling terms 

+ spin independent lead scattering-terms 

The first term resembles the Kondo model. The 
spin of the quantum dot Si = Y^aa 1 a L °Va-' a ia> 
couples to the band spin s^q 



The coupling strength 
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term describes a direct tunnel- 
quantum-dot electron to another 
effective coupling strength 



Ak 



+ 



+ 



with 



a denoting the opposite spin orientation of a. This term 
can also lead to a spin exchange in fourth order, so it 
is not negligible. The further parts of Eq. (|A2[) include 
processes, which change the occupation of one quantum 
dot by two electrons, and spin-independent scattering of 
continuum electrons at one dot. 



For the second Schrieffer- Wolff Transformation, the 
procedure is very similar. Using the ansatz for iSi that 
resembles the second order part of the interaction term, 
derived by the first transformation. However, as one is in 
the end interested in the fourth-order parts of the Hamil- 
tonian, which couple two quantum dot spins and conserve 
the quantum dot occupation number, one only needs to 
consider the first two parts of Eq. (|A2[) . The genera- 
tor for the transformation therefore can be written as 
1S2 



1S2 with 



^ = E 77 1 



i,kq,cr 



k — £q 



'JkqSi ' S kq 



(Bl) 



iS. 



(M 



- V 



e i + U i n«r-e j -U j n ji ," v l ° Ja 



The other second order terms finally drop out in the end, 
when the Hamiltonian is projected on the subspace of sin- 
gle occupied quantum dots. For the regrouping of terms 
in Eq. ([5][5]), one needs the symmetry of the expressions 
under the replacement 1 •<->• 2 and k -B- q. 
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